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A BLR origin for the iron Ka line in NGC 7213 
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ABSTRACT 

The X-ray spectrum of NGC 7213 is known to present no evidence for Compton 
reflection, a unique result among bright Seyfert Is. The observed neutral iron Kq line, 
therefore, cannot be associated with a Compton-thick material, like the disc or the 
torus, but is due to Compton-thin gas, with the Broad Line Region (BLR) as the most 
likely candidate. To check this hypothesis, a long Chandra HETG observation, together 
with a quasi-simultaneous optical spectroscopic observation at the ESO NTT EMMI 
were performed. We found that the iron line is resolved with a FWHM=2 400lgoo° 
km s _1 , in perfect agreement with the value measured for the broad component of 
the Ha, 2640 j^ km s" 1 . Therefore, NGC 7213 is the only Seyfert 1 galaxy whose 
iron Ko; line is unambiguously produced in the BLR. We also confirmed the presence 
of two ionised iron lines and studied them in greater detail than before. The resonant 
line is the dominant component in the Fe xxv triplet, therefore suggesting an origin 
in collisionally ionised gas. If this is the case, the blueshift of around 1 000 km s _1 of 
the two ionised iron lines could be the first measure of the velocity of a starburst wind 
from its X-ray emission. 
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1 INTRODUCTION 

A narrow component of the iron Kq line is almost in- 
variably present in Chandra high energy gratings and 
XMM-Newton CCD spectra of Active Galactic Nu clei (e.g. 
iBianchi et al, I l2007l : lYaaoob fe Padmanabhanll2004 and ref- 
erences therein). The line, typically unresolved with upper 
limits of several thousands km s _1 for the Full Width at Half 
Maximum (FWHM), must be produced far from the nu- 
cleus, either in th e torus envisaged in the Unification Model 
|Antonuccilll993l ) or in the Broad Line Region (BLR) . In the 
former case, and if the torus is Compton-thick, a Compton 
reflection component should also be present. In the latter 
much fainter Compton reflection component is ex- 
pected, and the intrinsic width of the iron line should be the 
same as that of the optical broad lines. The torus hypoth- 
esis seems to be pr e ferred for the majority of sources (e.g. 
IBianchi et all 12004 iNandral 120061 ) . but to unambiguously 
distinguish between these hypotheses high energy resolution 
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observations are required, which at present only Chandra 
can provide. 



The Seyfert 1/LINER NGC 7213 (z=0. 005839) presents 
a neg ligible amount of C o mpton reflection (R=Af2/27r < 
0.19: iBianchi et all 120031 . 12004 ) . This result, confirmed 
by three other BeppoSAX observations, is unique among 
bright Seyfert Is observed b y BeppoSAX l|Perola et al.ll2002l ; 
Risaliti 20021: lDadinall2008T ) . Therefore, the observed neutral 
iron line, whose equivalent width is ~ 80 eV, cannot be as- 
sociated with a Compton-thick material, like the disc or the 
torus, but is due to Compton-thin gas, like the BLR. The 
line width is unresolved by XMM-Newton EPIC pn, giving 
only an upper limit of about 8 000 km s _1 . In this Letter, we 
present a long Chandra grating observation of the iron line 
width, which turns out to be resolved and fully consistent 
with the FWHM of the broad component of the Ha line, as 
measured in a quasi-simultaneous optical observation. 
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2 DATA REDUCTION 
2.1 X-ray observation 

NGC 7213 was observed in two consecutiv e Chandra Ad- 
vance d CCD Imaging Spectrometer (ACIS: ICarmire et al.l 
2003) observations, with the H igh Energy Transmiss ion 
Grating Spectrometer (HETGS: ICanizares et al.l l2005h in 
the focal plane. The first observation (obsid 7742) was per- 
formed on 2007, August 6th, for 115 ks, while the second 
one (obsid 8590) on 2007, August 9th, for 35 ks. Data were 
reduced with t he Chandra Interactiv e Analysis of Obser- 
vations fCIAO: iFruscione et all l2006h 4.0.1 and the Chan- 
dra Calibration Data Base (CALDB) 3.4.2 software, adopt- 
ing standard procedures. First order High Energy Grating 
(HEG) and Medium Energy Grating (MEG) spectra were 
extracted for the source and the background. After having 
verified that there is no significant variability neither within 
any single observation nor between the two, we co-added 
them with the tool add_grating_Spectra, in order to have 
a single HEG spectrum and a single MEG spectrum, for a 

total exposure time of 148 ks each. 

All the fits were performed with Xspec 12.4.0 (jArnaudl 

1996), taking into account the instrumental resolution with 
the response matrix computed with standard procedures. 
Local fits were performed aroun d the i ron K g band with 
the unbinned spectra, using the ICashl ( 19761 ) statistic, in 
order to take advantage of the high spectral resolution of 
the gratings. Broad band (0.4-10 keV) fits were instead per- 
formed on the binned spectra, where the large number of 
counts in each bin (250) allowed the use of the \ 2 statistic. 
We did not analyse the 0th order data in this paper, since 
the count-rate per frame (~ 0.5 s _1 in the inner 2-pixel re- 
gion) is larger than the one suggested forpileup-free spectra 
(see The Chandra ABC Guide to Pileup^j, as confirmed by 
the resulting very flat spectrum. 



2.2 Optical observation 

In order to measure the physical width of the Ha line, we 
asked for director discretionary time (ID: 279.B-5054A) to 
obtain a quasi-simultaneous observation of NGC 7213 at the 
ESO NTT telescope (La Silla, Chile) . The observation was 
carried out in service mode on 2007, September 22th with 
the EMMI spectrograph equipped with grating #7 and a 1- 
arcsec width slit. The wavelength resolution was 2.49 A (the 
pixel scale was 0.82 A) in the range 5650-7140 A. Four spec- 
tra, with 30 sec long integrations each, were obtained with 15 
arcsec offsetted positions (pattern A B A B). The reduction 
process used standard MIDAS and IRAF facilities. The raw 
data were bias-subtracted, corrected for pixel-to-pixel varia- 
tions (flat field) and eventually sky-subtracted. Wavelength 
calibrations were carried out by comparison with exposures 
of He and Ar lamps, with an accuracy of 0.17 A (1 a). Rel- 
ative flux calibration was carried out by obser vations of the 
spectropho tometric standard star LTT9239 (|Hamuv et al.l 
Il992l . Il994h . The final S/N ratio was about 22. 
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3 SPECTRAL ANALYSIS 

In the following, errors correspond to the 90% confidence 
level for one interesting parameter (A^ 2 = 2.71), while er- 
ror bars in the plots correspond to 1 a. Energies and wave- 
lengths are reported in the rest-frame of the source. The 
Galactic column density along t he line of sight to NGC 72 13 
is included (2.04 x 10 20 cm" 2 : iDickev fc Lockmaij|l99ch . 



3.1 The Fe complex 

Figure [1] shows the HEG spectrum in the 6.-7.5 keV band. 
Three significant emission lines are apparent. The strongest 
is the neutral Ka line at 6.39 71° 'mi keV, consistent with 
iron less ionised than Fe XII (|House|[l96a i. The observed 

flux (2.9t° 9 7 P h cm " 2 s_1 ) and the EW (!20t3o eV ) are in 
agreement, within er rors, with those found by XMM-Newton 
|Bianchi et alj feOQS ) . The iron Ka line is actually composed 
of a doublet, K qi at 6.404 keV and Ka2 at 6.391 keV, with 
a flux ratio 2:1 (Bearden 1967). Although the separation of 
this doublet cannot be resolved by the Chandra HEG, we 
adopted a rigorous approach fitting two Gaussians, sepa- 
rated by 13 eV, the flux ratio fixed to the expected one and 
the width free to vary, but forced to be the same for the 
two lines. The width of the lines is resolved: a — 22lg° eV. 
Forcing the lines to be narrow [a = 0) leaves large residu- 
als and a significantly worse fit (AC = +15). We note here 
that the measure of the width is completely unaffected by 
this modelli ng with respec t to a single Gaussian, as already 
observed bv lYaqoob et al] |200ll ) for other sources. 

Let us investigate the origin of the observed broadening. 
A possible origin for the observed width may be the presence 
of a Compton Shoulder (CS). If the CS is modelled with a 
Ga ussian line w ith centroid energy at 6.3 keV and a = 40 
eV (|Mattll2002l 'l. while the iron line width is fixed to 0, the 
fit is significantly worse (AC = +14) and only an upper 
limit on the CS flux is found. This is not surprising, given 
the absence of a Compton reflection component associated 
to the CS. Another explanation for the width of the iron line 
can be blending with iron lines from higher ionisation states, 
but this would imply a larger value for the centroid energy of 
the resulting line and a clear asymmetric profile. Therefore, 
Doppler broadening is left as the most likely explanation 
for the width of the iron line, which would correspond to a 
FWHM=2400^ 1 00 km s" 1 (see Fig.©. 

Along with the neutral iron line, two other signifi- 
cant emission lines are observed in the HEG spectrum of 
NGC 7213, at 6.72ll°;°J° and 6.987±g;^g keV (see Fig. □ 
and Table[l|. The latter can be interpreted as emission from 
Fe xxvi, but its energy is blueshifted by a velocity of 900l|og 
km s" 1 with respect to the theoretical value of 6.966 keV. 
If the other emission line is identified with the forbidden 
component of Fe xxv Ka (6.637 keV), the blueshift must 
be 3 800lgQQ km s _1 . If instead it is identified with the res- 
onant component (6.700 keV), its corresponding blueshift 
would be 9OOI7QQ km s -1 , i.e. consistent with the one found 
for the Fe xxvi line. The fluxes and EWs of the two lines (see 
Table [p are consist ent with those found by XMM-Newton 
i|Bianchi et al]|2003l ). Finally, no neutral K/3 is found. The 
90% confidence level upper limit on its flux (4 x 10" 6 ph 
cm -2 s _1 ) leads to an upper limit on the ratio against Ka 
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Figure 1. Chandra HEG spectrum of NGC 7213 around the iron 
line, along with the best fit with three Gaussian emission lines. 
The spectrum is re-binned for clarity and the positions of the 
most important transitions are marked. 



Table 1. NGC 7213: the iron line complex in the Chandra HEG 
spectrum. 

Line E T AE FWHM Flux EW 

(1) (2) (3) (4) (5) (6) 



Fe i-xo 


6.400* 


q + 6 
3 -ll 


2400 +ii00 


2 9+°' 9 


120+ 40 


Fe xxv (r) 


6.700 


+21+ 10 


< 5500 


0.7±0.5 


24 ± 17 


Fe xxvi 


6.966* 


+21+ 19 


< 2700 


1.3 ±0.6 


60 ±30 



Notes - Col. (1) Identification . Col. (2) Theoretical energy in keV 
(NIST: iRalchenko et al]|2004h . *Weighted mean of the doublet. 
Col. (3) E obs -E T in eV, after correction for z=0.005839. Col. (4) 
km s" 1 Col. (5) 10~ 5 ph cm" 2 s" 1 . Col. (6) eV. 

Table 2. NGC 7213: optical emission lines in the ESO NTT 
EMMI spectrum. 

Line Aj AA FWHM Flux 

(1) (2) (3) (4) (5) 
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2*10- s 2.5*1CH 3x10- s 3.5*10- s 4x10- s 4.5x10 s 
Flux (ph cm- 2 s- 1 ) 



Figure 2. Flux vs. FWHM contour plot for the neutral iron line 
in NGC 7213. The curves refer to AC = 2.30, 4.61 and 9.21, i.e. 
confidence levels of 68, 90 and 99% for two interesting parameters. 

of 0.18, which is fully consistent with the expected value for 
neutral iron (see lMolendi et alJlioolh . 

While a broadband X-ray analysis of the Chandra spec- 
trum is beyond the scopes of this Letter, we note here 
that the 2-10 keV binned spectrum can be modelled by 
a simple powerlaw with V = 1.69 ± 0.01 and a flux of 
2.32 ± 0.04 x 10~ n erg cm -2 s , corresponding to a lu- 
minosity of ~ 1.7 x 10 42 erg s _1 . These values are perfectly 
consi stent with what fou nd by XMM-Newton 6 years be- 
fore l|Bianchi et al.l [200511. No other significant emission or 
absorption lines are found. 

3.2 The optical spectrum 

The optical spectrum of NGC 7213 was analysed between 
6300 and 6800 A. The most important atomic transitions 



[Oi] 


6300.32 


-0.3 


620+^ 


3.38 ±0.01 


[Oi] 


6363.81 


±2.8 


620+™ 


0.94 ±0.01 


Ha 


6562.79 


-1.1 


280+3°, 


-, no+0.13 
1 -°''-0.15 


Ho 


6562.79 


±2.6 


26401^ 


14 75+ 18 


[Nil] 


6548.06 


-0.1 


370 ± 30 


1.15 ±0.07 


[Nn] 


6583.39 


-0.1 


370 ± 30 


3.45 ±0.21 


[S ii] 


6716.42 


-0.2 


3501?° 


1 7QO+0-012 
x- ,oo _0.016 


[Sn] 


6730.78 


-0.3 


350+.?° 


i 77 o+0.011 
'°-0.016 



Notes.- C ol. (1) Identi fication. Col (2) Laboratory wavelength 
(A), in air (Bowcn 1960). Col. (3) \ h s ~ Mt after correction for 
z=0.005839. Typical error ~ ±0.5 A. Col. (4) km s" 1 (instru- 
mental resolution not removed). Col. (5) 10 — 14 erg cm -2 s — 1 . 

expected in this wavelength range are clearly observed in the 
spectrum (see Fig. [3}, so no attempt was made to remove the 
(unknown) contribution from the st ellar continuum. Note 
that t he optical spectrum analysed bv lFilippenko fe Halperr] 
(1984), to whom we will refer, was not galaxy-subtracted 
longward of 5500 A, either. 

Fig. [4] shows the Ha blended with the [N n] doublet. 
In order to disentangle any broad component of the Ha, 
we assumed that (i) F([Nll] A6583/F([N n] A6548)=3, as re- 
quired by the ratio of the respective Einstein coefficients, (ii) 
A2/A1 = 6583.39/6548.06 and (iii) the [Nil] lines are Gaus- 
sians with the same width. With these physical assumptions, 
a fit with three narrow emission lines (Ha plus the [N 11] 
doublet) is very poor and strongly suggests the presence of 
a broad component. The addition of this further compo- 
nent leads to a very good fit: see Tableland Fig. [3] for the 
4 Gaussian lines included in the fit and their convolution 
which nicely matches the observed profile. 

The widths of the narrow Ha and the [Nil] doublet 
are similar (see Table [2J|, while the FWHM of the broad 
component of Ha (2640lgQ° km s _1 ) clearly requires a dif- 
ferent origin, i.e. the BLR. Our results are qualitatively in 
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NGC7213: ESO IMTT/EMMI 




60 



S300 6400 6500 S600 6700 6800 

X (A) 

Figure 3. ESO NTT EMMI spectrum in the 6300 and 6800 A. The best fit is superimposed: the included lines are listed in Table [2] 



NGC7213: ESO NTT/EM Ml 




Figure 4. ESO NTT EMMI spectrum of the Ha- [Nil] complex. 
The best fit is the result of the convolution of the four Gaussian 
lines plotted with broken lines and listed in Table [2] 



agreement with those o f lFilippenko fc Hak>ernll|l984l ). in the 
sense that at least a broad component for the Ha is required 
by the data (they actually fit a total of seven components, 
but with the warning that the physical significance of all 
these lines is questionable). Moreover, the FWHM we de- 
rive for the Ha is close to the width of the broad compo nent 
of the H / 3 as m easured bv lFilippenko fc Halperr] (|l984T ) and 
IWinklerl (jl992h , around 4 000 and 3 200 km s" 1 , respectively. 

Finally, the [S n] doublet has a FWHM consistent with 
that of the [Nil] doublet and the narrow Ha component, 
while the [O i] doublet is broader, but all in agreement with 
an origin in the Narrow Line Region. 



4 DISCUSSION 

Differently from all the other Seyfert Is with broadband X- 
ray observations, NGC 7213 lacks a Compton reflection com- 
ponent, the clearest signature for the presence of a Compton- 
thick material. Therefore, the observed neutral iron line 
must originate in a Compton-thin material, be it the BLR 
or a Compton-thin torus. 

The Chandra HEG spectrum of NGC 7213 allowed us 
to resolve, for the first time, the iron line width at 2 400lg 1 g 
km s" 1 (FWHM). The analysis of a quasi-simultaneous opti- 



cal observation confirmed the presence of a broad component 
of the Ha line, for which we measured a FWHM=2640±go° 
km s _1 . The widths of the two lines are in very good agree- 
ment, which suggests that they are likely to be produced in 
the same material. 

To test if this scenario is possible, we have to verify if the 
observed EW of the iron line is in agreement with an origin 
in the BLR. Followi ng the detaile d procedure described by 
lYaaoob et ail l|200llL derived from lKrolik fc KallmarJ (|l987T ) 
adopting updated atomic data, the expected EW for the 
BLR in NGC 7213 can be written as: 



EW Fo ; 



34 



(0.35) ( 



10 23 cm- 



eV 



(1) 



This estimate assumes a spherically symmetric cloud 
distribution for the BLR, with a covering factor of f c , a col- 
umn density of A^h for each cloud and an iron abundance rel- 
ative to hydrogen of Af c = 4.68 x 10 -5 ll Anders fc Gre vessel 
1989). The powerlaw index of T = 1.69 derived from the 
Chandra and the XMM-Newton observations was further 
adopted in the derivation of |TJ. Assuming f c = 0.35, a 
column density around 3 x 10 23 cm -2 can reproduce an 
EW~ 100 eV, which is the order of magnitude found by 
Chandra and XMM-Newton. These values for f c and Nn are 
within th e ranges usua lly assumed in photoionization models 
of BLRs (lNetzerlll990ri. In particul ar, f c — 0.35 is the value 
derived bv lGoad fc Koratkarl i|l99Sl ) from observational con- 
straints in the Seyfert 1 NGC 5548, even if more 'canonical' 
values around 0.1 and 0.25 are generally found. The large 
covering factor needed to account for the iron line EW may 
be at odds with the regular intensity of the observed broad 
optical lines, bu t the weakness or ab sence of an UV bump in 
NGC 7213 (e.g. IStarling et~aT1l2005r i may result in a deficit 
of UV photoionizing photons and compensate the geomet- 
rical factor. In any case, / c and Nu can be lower, provided 
that is larger than solar and /o r the X-ray illumin ation 
of the BLR is anisotropic (see e.g. lYaaoob et al.l [200ll . and 
references therein). On the other hand, we cannot exclude a 
further contribution to the iron line EW from a Compton- 
thin torus, located on a pc-scale, but we stress here that 
we do not have any other piece of evidence to support its 
presence. 

The presence of emission lines from highly ionised iron 
in the Chandra HEG sp ectrum of NGC 721 3 confirms the 
XMM-Newton results bv lBianchi et al.l (120031 ). As suggested 
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by these authors, their origin may be in gas photoionised 
by the AGN, as found in many Seyfert Is and 2s (e.g. 
IBianchi fe Matd 120021 : iBianchi et aD [2005). However, the 
gratings data presented here revealed that the resonant line 
is the dominant component in the Fe xxv triplet (unless its 
blueshift is significantly larger than the one measured for 
the Fe xxvi line), which is suggestive of an origin in gas in 
collisiona l equilibrium. This is in agreement with the results 
found bv lStarling et all (|2005l ) using diagnostic tools based 
on the O VII triplet observed in the XMM- Newton RGS 
spectrum of NGC 7213. It is interesting to note that the 
Fe xxv-xxvi lines' blueshift of around 1000 km s -1 (con- 
sistent wi th the centroid ener gies of the O VII lines as mea- 
sured by IStarling et ail [2005). together with the evidence 
that the gas producing such lines is in collisional ionization 
equilibrium, suggest that the hot, line emi tting, gas be in 
the form of a starburst driven wind (see e.g. Heckmanlliool 
and references therein), rather than an intrinsic (typically 
photoionized) AGN wind. 

Th e Ha morphology of NGC 7213 was studied in de- 
tail b y IStorchi-Bergmann et al.l (| 19961 ) and lHameed et al.l 
(2001). A circumnuclear ring of star formation is located 
at 20-30 arcsec from the nucleus (corresponding to a 
few kpc) and a giant Ha filament is present at 2.9 ar- 
cmin of distance (around 19 kpc). Neither of these struc- 
tures is observed in the Chandra X-ray image, which only 
shows an unresolved (wi t hin fe w arcsec) nucleus. However, 
IStorchi-Bergmann et al.l (| 19961 ) reported the possible pres- 
ence of a collimated outflow within 14 arcsec from the nu- 
cleus, with a velocity around 100 km s _1 . The gas producing 
the ionised iron lines reported in this Letter may be the in- 
ner, hotter and faster phase of this superwind. If this is the 
case, this could be the first measure of the velocity of a 
starburst wind from its X-ray emission. 
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